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INTRODUCTION 


The  electroencephalograph  (EEG)  consists  of  an  ensemble  of  electric  potentials  recorded 
continuously  from  a  standard  montage  of  scalp  electrodes  and  is  a  direct  monitor  of  the  central 
nervous  system  (CNS).  During  surgery  the  EEG  is  used  to  correlate  the  activity  of  the  CNS  to  depth 
of  anesthesia.  Depth  of  anesthesia  is  related  to  the  patient’s  ability  to  perceive  and  respond  to  noxious 
stimuli.  Historically,  anesthesiologists  have  used  a  variety  of  examination  procedures  to  estimate 
awareness  of  a  subject  (ref  1).  The  use  of  electroencephalography  to  monitor  depth  of  anesthesia  was 
first  suggested  by  Gibbs  et  al.  (ref  2),  who  observed  the  correlation  between  the  EEG  patterns  and 
levels  of  anesthesia  induced  by  diethyl  ether.  With  the  development  of  fast  and  inexpensive 
computers,  various  signal  processing  techniques  have  been  utilized  to  analyze  and  quantify  the  signals 
both  in  the  time  or  in  the  frequency  domaia  Klein  and  Davis  (ref  3)  used  the  frequency  and 
amplitude  of  the  zeroth,  first,  and  second  derivative  of  the  signal  to  monitor  patients  subjected  to 
general  anesthesia.  A  hierarchical  clustering  identification  system  based  on  autoregressive  (AR) 
estimation  of  the  EEG  spectra  was  designed  by  Thomsen  et  al.  (ref  4).  Autoregressive  estimation  was 
used  by  Cerruti  et  al.  (ref  5)  in  monitoring  the  influence  of  hypotension  on  the  CNS  during  riskful 
neurosurgery.  A  tenth  order  AR  model  was  also  utilized  by  Sharma  and  Roy  (ref  6)  to  train  a  neural 
network  to  predict  the  anesthetic  depth  (ref  6). 

In  the  frequency  domain,  Fourier  analysis  has  been  used  to  measure  the  spectral  properties  of 
the  EEG.  The  spectral  information  (ref  7)  is  usually  obtained  in  terms  of  the  median  frequency, 
special  edge  frequency  (SEF),  total  power,  and  relative  power  in  the  individual  frequency  bands. 

These  bands  have  been  defined  as  Delta  (1.0  to  3.5  Hz),  Theta  (3.5  to  7.5  Hz),  Alpha  (7.5  to  12.5  Hz), 
Betal  (12.5  to  17.5  Hz),  and  Beta2  (17.5  to  25  Hz).  In  the  Fourier  transform,  the  signal  is  assumed  to 
be  stationary  in  the  processing  period,  i.e.,  it  is  assumed  that  its  spectrum  does  not  change  over  time. 

It  was  shown  by  Koch  et  al.  (ref  8)  that  the  EEG  spectral  parameter  SEF  can  lead  to  contradictory 
results,  since  it  is  based  on  the  assumption  of  stationarity.  The  authors  observed  that  the  SEF  not  only 
decreased  when  the  patient  was  deeply  anesthetized,  but  also  when  surgically  stimulated  while  under 
light  anesthesia.  Such  occurrences  have  been  labeled  as  "paradoxical"  arousal.  For  this  reason,  time- 
frequency  techniques  have  been  developed  to  study  the  spectral  changes  in  the  EEG.  The  most 
common  technique  is  the  Short-Time  Fourier  Transform  (STFT)  (ref  9).  It  is  obtained  by  sliding  a 
window  across  the  signal  and  computing  the  magnitude  of  the  Fourier  transform  of  the  windowed 
segment  with  the  assumption  that  the  signal  is  stationary  in  the  time  interval  defined  by  the  window. 
The  spectral  resolution  of  the  STFT  can  be  improved  by  using  a  longer  window.  A  longer  window, 
however,  leads  to  smearing  of  nonstationary  data,  while  a  shorter  window  leads  to  poor  frequency 
resolution.  Therefore,  Nayak  et  al.  (ref  10)  studied  the  effect  of  noxious  stimuli  on  the  EEG  by  using 
a  time-frequency  spectral  representation  (TESR)  known  as  the  Choi-Williams  Exponential  Distribution 
(ED)  (ref  1 1).  The  TFSRs  during  stimuli-induced  wakefulness  were  localized  in  time  and  frequency 
with  increased  Delta  activity.  In  deeply  anesthetized  states,  although  the  frequencies  were  in  fte  Delta 
band,  the  spectra  were  not  localized.  Although  the  results  gave  a  broad  picture  of  changes  in  the 
EEG,  finer  details  could  not  be  obtained  due  to  the  presence  of  crossterms  (ref  12).  These  crossterms 
carry  redundant  information  and  obscure  important  features  of  the  signal. 

The  wavelet  transform  (WT)  is  the  latest  technique  for  processing  signals  with  nonstationary 
spectra  (refs  13,14).  Applications  range  in  many  fields,  such  as  ultrasonics  (ref  15),  geophysics, 
mathematics,  theoretical  physics,  and  communication  (ref  16).  Multiresolution  signal  decomposition 
using  wavelet  transform  has  been  extensively  studied,  especially  in  the  development  of  Perfect 
Reconstruction-Quadrature  Mirror  Filter  (PR-QMF)  (ref  17). 
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The  WT  is  defined  in  tenns  of  basis  functions  obtained  by  compression/dilation  and  shifting  of 
a  "mother  wavelet."  It  does  not  suffer  from  the  time-bandwidth  resolution  tradeoff  like  the  STFT. 
Since  it  is  a  linear  transform,  the  crossterms  associated  with  the  ED  are  avoided  (ref  18).  These 
features  make  the  WT  an  attractive  tool  to  study  nonstationary  information  in  the  EEG.  The  ability  of 
the  WT  to  use  different  windows  at  different  frequencies  would  tend  to  model  the  EEG  in  an  efficient 
way.  The  windows  at  different  frequencies  have  a  constant  Q  factor  =  isflj),  where /is  the  center 
frequency  and  A/"  is  the  spread  of  the  window  spectrum.  Therefore,  the  stimuli-induced  transients  in 
the  EEG  can  be  studied  more  effectively  by  using  the  WT  and  lead  to  possible  qualification  of  the 
EEG  spectral  information  in  relation  to  anesthetic  depth. 

In  the  following  the  wavelet  transform  is  briefly  introduced  and  a  Wavelet  Transfonn  .Signal 
Processor  is  designed  for  the  investigation  of  the  EEG  signals.  The  choice  of  the  mother  wavelet 
defines  the  type  and  class  of  information  that  can  be  extracted.  In  this  work  the  EEG  signals  were 
analyzed  using  two  different  mother  wavelets.  In  one  case  the  temporal  resolution  of  the  EEG 
spectrum  is  studied,  while  the  other  mother  wavelet  is  used  in  detection  and  analysis  of  transient  and 
background  processes  in  the  brain.  The  advantages  and  disadvantages  of  each  choice  are  presented  in 
connection  with  the  information  obtained  in  each  case. 

SIGNAL  PROCESSING  USING  WAVELET  TRANSFORM 

By  definition  the  WT  is  the  correlation  between  the  signal  and  a  set  of  basic  wavelets  (ref  19). 
An  appropriate  square  integrable  mother  wavelet  hif)  is  chosen  to  analyze  a  specific  transient  signal  of 
finite  energy.  Then  a  complete  orthogonal  set  of  daughter  wavelets  is  generated  from  the 
mother  wavelet  h(t)  by  dUation  {a)  and  {b)  shift  operations.  The  WT  expansion  coefficient  WXa,b)  of 
the  signal  s{t)  are  given  by 

-Q* 

=  s(f)0— A’fl 

where  the  function  is  given  by 

Kt<f>  ’ 

The  normalization  constant  a  is  chosen  so  every  daughter  wavelet  has  identical  energy.  The  Fourier 
transform  of  the  daughter  wavelet  h^j,(t)  is  given  by 

=  v/5  •  (3) 

where  H(f)  represents  the  Fourier  transform  of  the  mother  wavelet  This  equation  shows  the  important 
concept  that  a  dilation  t/a  in  the  time  domain  is  equivalent  to  a  frequency  change  of  (rf,  while  the  time 
delay  b  results  in  a  phase  shift  in  the  frequency  domain  and  does  not  affect  the  magnitude  of  the 
Fourier  spectrum.  If  the  mother  wavelet  can  be  represented  as  a  pass-band  filter  of  central  frequency 
/,  and  bandwidth  A/,,  then  the  daughter  wavelet  is  a  pass-band  filter  of  a*/,  and  of  central 
frequency  and  bandwidth,  respectively.  The  relative  bandwidth //A/’=  Q  is  thus  constant  and  does  not 
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depend  on  the  value  of  the  dilation  a.  The  WT  can  thus  be  seen  as  a  bank  of  filters  constructed  by 
dilation/compression  of  a  single  function  h(t).  The  filter  bank  structure  and  the  relationship  between 
the  output  WXa,b)  and  the  two  inputs,  sO)  and  hit),  are  shown  clearly  in  Figure  1.  The  output  is  not 
only  related  to  the  input  sit),  but  also  to  the  type  of  filters  h^it)  that  are  utilized.  By  choosing  a 
different  mother  wavelet,  different  characteristics  of  flie  input  signal  sit)  can  be  obtained  as  output. 

The  filter  constructed  by  the  dilated  version  of  the  mother  wavelet  processes  the  low  frequency 
information  of  the  signal  sit),  and  the  one  related  to  the  compressed  version  of  hit)  analyzes  the  high 
frequency  information.  The  ability  of  adapting  the  window  size  to  the  signal  being  processed  makes 
the  WT  a  natural  candidate  for  analysis  of  the  transient  waveforms  with  a  wide  spectral  range  (ref  20), 

The  WT  has  been  introduced  as  an  improvement  over  the  windowed  or  STFT 
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where  wit)  is  an  appropriate  window  function  in  the  time  domain.  If  the  window  w(r)  is  Gaussian  in 
shape,  it  is  commonly  referred  to  as  the  Gabor  transforai  (ref  21).  The  STFT  analysis  has  the  same 
resolution  at  all  locations  in  the  time-frequency  plane  due  to  the  fact  that  the  same  time  window,  wit), 
is  used  over  the  entire  frequency  range.  In  the  WT  the  frequency  resolution,  and  thus  also  the 
temporal,  is  frequency-dependent,  resulting  in  a  finer  frequency  resolution  in  the  lower  range  and  a 
finer  temporal  resolution  in  the  higher  frequency  range.  This  property  of  the  WT,  referred  to  as  the 
constant  Q  factor,  results  in  a  more  dynamic  representation  of  the  signal  in  the  time-frequency  or  time- 
scale  plane. 

Another  difference  from  the  STFT  is  that  the  mother  wavelet  can  be  any  function  that  satisfies 
the  following  admissibility  condition  (ref  22): 

rfco  =  c  <  -  (s) 

0  " 

where  co  =  lirf  is  the  angular  frequency. 

In  the  Fourier  analysis  the  signal  is  expanded  into  a  series  of  complex  exponentials,  while  in 
the  WT  the  signal  is  decomposed  into  a  series  of  elementary  wavelets,  which  can  be  modeled  to  match 
the  response  of  a  system  to  external  disturbances.  The  flexibility  of  choosing  the  proper  mother 
wavelet  is  one  of  the  strongest  advantages  of  using  the  WT,  since  the  choice  of  mother  wavelet  for  a 
particular  problem  improves  the  signal  processing  capability  of  the  technique.  Tailoring  of  the  wavelet 
is  possible  and  should  be  done. 

The  Morlet  and  the  Mexican  hat  mother  wavelets  are  popular  in  the  signal  processing 
community.  The  Morlet  wavelet  is  expressed  as 
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where  A  is  the  amplitude,  x  is  the  decay  constant,  ©  is  the  angular  frequency,  and  (j)  is  the  phase  shift. 
The  real  and  complex  parts  of  the  Morlet  wavelet  are  shown  in  Figure  2a  and  its  spectra  in  Figure  2b. 
It  can  be  seen  that  the  spectra  of  the  wavelet  can  be  represented  by  a  Gaussian  distribution  with  center 
frequency  at  0.5  Hz  and  a  j2  equal  to  5.  The  mother  wavelet  in  Eq.  (5)  is  represented  by  the  product 
of  two  components:  a  high  frequency  complex  harmonic  component  and  a  low  frequency  real 
Gaussian  pulse.  The  high  ftequency  component  is  intentionally  constructed  to  be  a  complex  function. 
This  eliminates  oscillations  in  the  plot  of  magnitude  versus  time.  The  Gaussian  nature  of  the  function 
in  time  as  well  as  frequency  domain  has  proved  to  be  useful  in  studying  the  spectral  evolution  of 
signals  over  time. 

The  Mexican  hat  wavelet  is  used  in  analyzing  pulses  in  signals.  Pulse  detection  consists  of 
determining  the  presence  or  absence  of  a  pulse  and  estimating  its  amplitude  and  arrival  time  even 
when  the  pulse  is  buried  in  time.  The  real  part  of  the  Mexican  hat  function  is  given  by 


=  A 

1  - 

f-Tl 

j 

(7) 


where  A,  u,  and  v  are  real  constants.  Usually,  a  complex  Mexican  hat  mother  wavelet  is  used  to 
suppress  oscillations  in  the  magnitude  plot.  When  detection  of  pulse  is  of  primary  concern,  the 
parameters  A,  u,  and  v  are  chosen  so  as  to  have  a  Q  factor  equal  to  one.  The  real  and  imaginary  parts 
of  such  a  Mexican  hat  mother  wavelet  and  its  spectrum  are  shown  in  Figure  3.  As  the  Q  factor  is 
equal  to  one,  the  processes  are  depicted  in  different  dilation  ranges  as  long  and  short  responses. 
Therefore,  the  Mexican  hat  approach  with  Q  equal  to  one  can  be  used  to  analyze  the  signal  structure 
in  the  time  domain. 

In  many  practical  applications,  as  the  one  discussed  in  this  report,  the  signal  is  digitized  in 
discrete  time  intervals.  Furthermore,  as  for  the  STFT,  the  WT  is  highly  redundant  when  the 
parameters  (i,b  are  continuous  in  value.  Therefore,  the  WT  is  evaluated  for  a  discrete  number  of 
points  in  the  time-scale  plane,  corresponding  to  a  discrete  set  of  wavelets.  Discretization  of  the  time- 
scale  parameters  (ref  23)  leads  also  to  setting  tz  =  2"  and  b  =  n  ’  2",  where  it  is  assumed  unity  for  the 
discrete  sampling  interval.  The  discrete  WT  is  thus  expressed  as 

<*) 

*=i 


DATA  COLLECTION  AND  DATA  PROCESSING 

The  EEG  signals  shown  in  Figure  4  were  obtained  from  mongrel  conditioned  dogs  weighing 
approximately  20  to  25  Kg.  The  dogs  were  anesthetized  by  orally  inmbating  2  mg/Kg  doses  of 
Methohexital  (Brevital).  The  subjects  were  placed  on  a  closed-circuit  anesthesia  system  (ref  24)  at  an 
anesthetic  level  high  enough  for  the  placement  of  monitors  without  animal  movement. 
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Electrocardiograph  leads  were  placed  using  needle  electrodes,  and  four  channels  of  EEG  needle 
electrodes  (ref  25)  were  applied  as  shown  in  Figure  5.  The  four  channels  were  LfL^,  Rf-Ro>  Fz-Oz>  V- 
R,.  Arterial  pressure  was  measured  using  a  percutaneously  placed  femoral  arterial  line.  Arterial 
pressure  and  electrocardiogram  (ECG)  signals  were  analyzed  using  a  Mermen  Medical  Horizon 
monitor,  and  the  EEG  signal  was  monitored  by  an  Axon  System  Sentinel-4EEG/EP  monitor. 
Respiratory  gases  were  monitored  by  a  PPG  Medspect  mass  spectrometer  and  a  Criticare  Systems 
POET  capnograph.  The  capnograph  provided  a  continuous  visualization  of  the  subjects’  breathing 
pattern.  All  signals  were  processed  by  a  486  PC  computer  equipped  with  a  Data  Translation  DT2801 
analog-to-digital  acquisition  board. 

The  anesthetic  agent  (halothane)  was  increased  fiom  0.2  to  1.2  percent  in  increments  of  0.2 
percent  at  time  intervals  of  approximately  30  minutes.  At  the  end  of  each  interval,  when  the 
anesthetic  level  had  stabilized,  as  determined  by  the  mass  spectrometer,  the  subjects  were  tested  for 
depth  of  anesthesia.  This  test  was  accomplished  by  tail  clamping  for  30  seconds.  This  results  in  a 
supramaximal  stimulus  (ref  26).  At  the  time  of  tail  clamping,  the  depth  of  anesthesia  was  estimated 
by  observation  of  animal  movement  and  changes  in  respiration,  heart  rate,  and  mean  arterial  pressure. 
The  depth  was  graded  as  either  awake  or  asleep.  At  the  conclusion  of  the  experiment,  the  animal  was 
awakened,  extubated,  and  returned  to  the  animal  facility.  This  protocol  was  approved  by  the  Institute 
of  Animal  Care  and  Use  Committee. 

A  Wavelet  Transform  Signal  Processor  was  designed  to  investigate  the  EEG  signals.  The 
system,  shown  in  Figure  1,  receives  the  EEG  signal  as  input,  h(t)  is  the  mother  wavelet,  and  W,(a,b)  is 
the  WT  of  EEG.  The  WT  technique  was  applied  to  the  EEG  signals  obtained  ftom  Channel  I.  Two 
different  mother  wavelets  were  used  for  analysis  of  the  EEG.  In  one  case,  the  Morlet  function  was 
used  as  the  mother  wavelet  to  optimize  the  resolution  of  time-scale  analysis.  A  Mexican  hat  function 
was  obtained  by  fitting  the  EEG  response  for  the  awake  subject  during  stimulus  and  was  used  as  the 
mother  wavelet  to  detect  the  presence  of  spikes  corresponding  to  the  EEG  obtained  during  light 
anesthesia.  Thus,  the  WT  transform  can  be  interpreted  as  the  correlation  or  affinity  between  the  EEG 
signal  and  the  brain  response  to  stimulus  modeled  by  the  Mexican  hat  shown  in  Figure  3.  The  EEG 
data  before  and  during  the  period  of  tail  clamping  were  analyzed. 

TIME-SCALE  ANALYSIS  OF  EEG  USING  WAVELET  TRANSFORM 

The  spectral  properties  of  the  EEG  signals  were  analyzed  in  the  frequency  range  from  d.c.  to 
25  Hz.  The  center  frequency/,  of  h,(t)  is  0.5  Hz,  therefore  for  a  generic  dilation  m,  the  center 
frequency  of  the  associated  filter  is  0.5/2'".  In  order  to  cover  the  required  range  of  frequencies,  the 
coefficient  m  has  to  be  as  low  as  -6(f„  =  0.5/2"*  =  32  Hz)  and  as  high  as  +l(/„  =  0.5/2*  =  0.25  Hz). 
Therefore,  the  absolute  bandwidth  will  vary  from  6  to  0.05  Hz  for  m  =  -6  and  m  =  1,  respectively. 

The  time-scale  representation  of  the  WT  for  the  EEG  obtained  during  state  of  light  (depth  zero)  and 
deep  anesthesia  (depth  one)  are  given  as  density  plots.  The  x  and  y-axes  in  the  plot  represent  the 
dilation  coefficient  m  and  the  time  coefficient  n,  respectively.  A  negative  dilation  coefficient  m 
represents  compression  of  the  mother  wavelet  and  thus  higher  frequencies.  The  correspondence 
between  dilation  coefficient  m  and  the  frequency  bands  of  interest  is  given  in  Table  1.  The  magnitude 
of  the  WT  is  color-coded  in  a  logarithmic  scale  in  order  to  show  contributions  to  the  EEG  signals 
from  the  various  frequency  bands.  Red  corresponds  to  higher  values  of  the  WT  coefficients  and  blue 
corresponds  to  lower  values.  Also  in  each  case,  depth  zero  and  one,  the  magnitude  of  the  WT 
coefficients  was  time-averaged  over  the  pre-stimulus  period  and  the  stimulus  period.  These  are  plotted 
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as  a  function  of  the  center  frequency  of  the  filter  associated  with  each  daughter  wavelet.  Thirteen 
signals  were  processed  using  the  WT-seven  were  classified  as  asleep  and  five  were  classified  as 
awake.  Typical  awake  and  asleep  signals  are  discussed  below. 

Figure  6  shows  the  density  of  plot  of  the  time-scale  representation  when  depth  was  equal  to 
zero.  A  local  maximum  of  the  WT  can  be  observed  before  the  stimulus  for  m  =  -5.5  (22.6  Hz),  which 
seems  to  shift  to  m  =  -5.75  (29.6  Hz)  once  the  stimulus  is  applied.  A  somewhat  lower  density  of 
values  is  observed  during  the  stimulus  in  the  range  m  =  -5  to  m  =  -3  (4  to  16  Hz),  but  it  is  difficult  to 
quantify  the  variations.  From  the  plot  of  the  averages  of  the  wavelet  coefficients.  Figure  7,  it  can  be 
seen  that  there  is  a  dramatic  increase  in  magnitude  of  the  low  frequency  WT  coefficients  during  the 
stimulus  period.  The  shift  in  the  Beta2  range,  observed  in  the  previous  plot,  is  also  detected.  A 
decrease  in  ampUtude  is  also  observed  in  the  Alpha  range. 

The  time-scale  representation  of  the  WT  for  the  sleep  signal,  depth  one,  is  shown  in  Figure  8 
as  a  density  plot.  The  wavelet  magnitudes  are  much  lower  than  those  given  in  Figure  6.  The  large 
changes  in  amplitude  observed  in  the  awake  case  are  not  present  here,  in  particular  the  lower 
frequency  activity  seems  to  be  unrelated  to  the  stimulus.  But  some  shifts  can  be  observed  in  the 
higher  frequency  ranges;  in  particular  the  shift  between  m  =  -3.7  to  m  =  -4.1  (6.5  to  8.6  Hz),  which 
corresponds  to  a  Theta-to-Alpha  region  shift,  and  m  =  -4.7  to  m  =  -5.1  (13  to  17  Hz)  in  the  Betal 
range.  These  observations  are  confirmed  by  the  one-dimensional  plot  of  Figure  9.  A  more  subtle  shift 
of  values  can  be  observed  in  the  m  =  -2  to  m  =  -3  range  (2  to  4  Hz),  which  is  the  Delta  region,  but  it 
is  not  possible  to  quantify  the  variations. 

SYSTEM  RESPONSE  ANALYSIS 

Figures  10  and  11  show  the  corresponding  WT  for  the  awake  and  asleep  signals,  respectively. 
As  expected,  the  WT  for  the  awake  signal  shows  that  the  Mexican  hat  is  a  strong  indicator  of 
stimulated  activity,  which  again  is  concentrated  in  the  high-scale  Gower  frequencies)  range. 
Interestingly,  for  the  sleep  signal,  the  WT  indicates  that  the  stimulus  depresses  the  system  response  for 
dilation  values  between  m  =  -\  to  /n  =  1. 

A  one-dimensional  plot  of  with  the  dilation  m  set  to  l(a  =  0.5)  is  shown  in  Figures 
12a  and  12b  for  awake  and  asleep  states.  It  can  be  observed  that  the  response  of  both  states  before  the 
appUcation  of  the  stimulus  is  similar  in  magnitude.  However,  the  magnitudes  are  higher  for  the  depth 
zero  EEG  when  compared  to  those  at  depth  one. 

DISCUSSION 

From  the  time-scale  analysis  of  the  EEG  using  the  Morlet  mother  wavelet,  it  can  be  concluded 
that  for  the  sleep  case,  a  shift  toward  higher  frequencies  can  be  observed  in  the  Betal  region  which  is 
analogous  to  the  shift  observed  in  the  Beta2  region  for  the  awake  case.  The  average  magnitudes  of 
the  pre-stimulus  WT  coefficients  corresponding  to  different  frequencies  for  the  unstimulated  sleep  and 
awake  states  are  plotted  in  Figure  13.  The  WT  coefficients  have  roughly  the  same  order  of  magnitude 
in  both  cases,  but  they  show  different  spectra,  especially  in  the  higher  frequency  range.  The  awake 
state  is  characterized  by  Beta2  (22.5  Hz)  waves  and  lower-amplitude  Alpha  waves  (7.5  Hz),  while  the 
sleep  state  is  characterized  by  Betal  (15  Hz)  waves  and  slightly  higher-amplitude  Alpha-Theta  (5  to  8 
Hz)  waves.  It  has  been  suggested  that  Alpha  blocking  can  be  related  to  the  state  of  the  patients  (ref 
27)  and  thus  can  be  utUized  to  evaluate  the  depth  of  anesthesia.  Alpha  blocking  is  a  response  to 
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noxious  stimulus  under  light  anesthesia  conditions.  Our  results  show  the  presence  of  Beta2  shift  in 
addition  to  Alpha  blocking.  During  deep  anesthesia  a  Betal  shift  seems  to  exist.  There  also  seems  to 
be  three  distinct  conditions  in  which  the  subject  can  be  found:  (1)  asleep  with  higher- amplitude  Alpha- 
Theta  and  Betal  waves;  (2)  awake,  with  lower-amplitude  Alpha-Theta  waves  and  significant  Beta2 
activity;  and  (3)  stimulated  awake  state,  with  a  large  increase  in  amplitude  in  the  low  frequency  range 
and  a  shift  in  the  main  fi-equency  in  the  Beta2  region. 

The  Beta2  shift  in  the  case  of  depth  zero  EEG  and  Betal  shift  associated  with  depth  one  EEG 
could  not  be  detected  by  the  system  analysis  technique  using  the  Mexican  hat  wavelet.  This  could 
suggest  that  the  Beta  fr^uency  shifts  are  not  due  to  external  stimulus  but  due  to  some  internal 
process.  The  difference  in  response  at  depth  zero  and  the  one  due  to  the  presence  of  stimulus  suggests 
a  background  activity  of  the  brain.  The  brain  reacts  to  the  external  stimulus  in  two  different  modes 
depending  on  the  state  of  consciousness  of  the  subject.  In  the  awake  case,  there  is  an  evident  increase 
in  response,  while  for  sleep  a  reduction  in  this  activity  is  observed.  This  analysis  seems  to  suggest 
that  the  brain  has  an  ongoing  background  process  that  monitors  external  stimulus  in  both  the  sleep  and 
awake  states.  Similar  conclusions  were  reached  by  authors  in  Reference  10. 

COMMENTS  AND  CONCLUSIONS 

The  wavelet  transform  was  introduced  and  used  for  the  spectral  analysis  of  the  EEG  signals  in 
the  awake  and  asleep  states.  The  choice  of  the  mother  wavelet  used  in  the  analysis  defines  the  type  of 
information  obtained  by  such  analysis,  for  this  reason  two  different  mother  wavelets  were  utilized  in 
this  work.  Using  the  WT  with  the  Morlet  wavelet,  the  awake  and  asleep  states  were  characterized  by 
the  frequency  contents.  This  approach  is  similar  to  the  STFT.  Differences  in  the  frequency  spectra 
between  the  two  states  were  observed  and  are  summarized  as  follows:  (1)  for  the  awake  state,  the 
stimulus  causes  an  increase  in  low  frequency  energy;  (2)  for  the  asleep  state,  stimulus  causes  the  dog’s 
EEG  waves  to  shift  in  the  awake  direction,  resulting  in  higher  frequency  Alpha  waves  and  a  shift  from 
Betal  to  Beta2  waves;  (3)  there  exists  a  possibility  of  predicting  the  state  of  a  subject  by  recording  the 
changes  in  Alpha  and  Beta  waves. 

From  these  observations,  there  seems  to  be  three  distinct  conditions  with  probably  some 
transition  states  in-between.  These  three  states  are  (1)  sleep,  with  the  spectrum  mainly  localized  in  the 
Alpha-Theta  and  Betal  regions;  (2)  awake,  with  the  spectrum  localized  in  the  lower  frequency  Alpha 
region  and  in  the  Beta2  range;  and  (3)  stimulated  awake  state,  a  dramatic  increase  in  low  frequency 
energy  on  top  of  the  regular  awake  frequency  spectrum. 

The  second  choice  of  the  mother  wavelet  was  to  model  the  brain  response  to  the  external 
stimulus.  Using  the  EEG  signal  in  the  awake  state,  a  Mexican  hat  mother  wavelet  was  constracted 
and  used  to  study  the  EEG  signals.  Results  show  that  there  seems  to  be  a  constant  background 
process  that  monitors  internal  and  external  stimuli.  This  background  process  is  the  same  whether  in 
the  awake  or  sleep  states  if  there  is  not  supramaximal  external  stimulus.  In  the  awake  state,  this 
external  stimulus  increases  the  system  response;  while  in  the  asleep  state,  the  opposite  seems  to 
happen,  resulting  in  a  disruption  of  the  background  process.  This  approach  modeled  in  Figure  14 
could  be  used  to  model  the  brain  functions  using  EEG  data. 
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Table  1.  Correspondence  Between  the  Frequency  Bands  and  the  Dilation 
Coefficients  m  Used  in  the  Wavelet  Transform  Analysis  of  the  EEC  Signals 


Band 

Frequency  (Hz) 

Dilation  m 

Delta  (A5) 

1.0-3.5 

1. 0-2.8 

Theta  (00) 

3.5-7.5 

2.8-3.9 

Alpha  (A.a) 

7.5-12.5 

3.9-4.64 

Betal  ((31) 

12.5-17.5 

4.64-5.13 

Beta2  (p2) 

17.5-25 

5.13-5.64 
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h(t) 

Figure  1.  Block  diagram  of  the  Wavelet  Transform  Signal  Processor  designed  for  the 

investigation  of  the  EEG  signals.  The  top  plot  shows  the  two  inputs  s(f)  and  h{t) 
and  the  relative  output  W,(a,b%  while  the  expanded  view  in  the  tottom  shows  the 
filter  bank  structure  of  the  processor  and  the  relationship  of  the  output  with  the  inputs. 
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Magnitude  Amplitude 


Magnitude  Amplitude 


Frequency  [Hz] 


Figure  3.  Time  and  frequency  plots  of  the  Mexican  hat  mother  wavelet  of  Eq.  (7)  used  to 
model  the  brain  response  to  an  external  stimulus.  The  real  and  complex  parts 
of  the  wavelet  are  plotted  in  Figure  3a,  and  its  spectra  are  given  in  Figure  3b. 
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Figure  4.  EEG  signals  obtained  from  the  subjects  in  the  awake  and  sleep  state, 
respectively.  The  square  signal,  overlapped  to  the  experimental  data, 
is  high  when  the  dog  was  subjected  to  the  stimulus. 
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Fz  -  Frontal  z 
Lf  -  Left  Frontal 
Rf  -  Riffht  Frontal 
Lo  -  Left  Occipital 
Ro  -  Right  Occipital 
Ox  -  Occipital  z 
Rt  -  Right  Temporal 
V  -  Vertex 
Gnd  -  Ground 


Figure  5.  Electrode  configuration  used  in  the  EEG  experiments. 
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WT  Awake  Signal  [Log  scale] 
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Figure  6.  Time-scale  representation  of  the  WT  for  the  awake  signal  using  the  Morlet 
mother  wavelet  hjit).  The  x  and  y-axes,  respectively,  represent  the 
dilation  coefficient  m  and  the  time  delay  coefficient  n,  as  defined  in 
Eg.  (5).  The  magnitude  of  the  WT  is  color-coded  in  a  logarithmic  scale. 
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WT  magnitude 


17 


WT  Sleep  Signal  [Log  scale] 
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Figure  8.  Time-scale  representation  of  the  WT  for  the  sleep  signal  using  the  Morlet 
mother  wavelet  /i/r).  The  x  and  y-axes,  respectively,  represent  the 
dilation  coefficient  m  and  the  time  delay  coefficient  n,  as  defined  in 
Eq.  (5).  The  magnitude  of  the  WT  is  color-coded  in  a  logarithmic  scale. 
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Figure  10.  Time-scale  representation  of  the  WT  for  the  awake  signal  using  the  Mexican 
hat  mother  wavelet  /tjCr).  The  x  and  y-axes,  respectively  represent  the 
dilation  coefficient  m  and  the  time  delay  coefficient  n,  as  defined  in 
Eq,  (5).  The  magnitude  of  the  WT  is  color-coded  in  a  logarithmic  scale. 
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Figure  11.  Time-scale  representation  of  the  WT  for  the  sleep  signal  using  the  Mexican 
hat  mother  wavelet  /ijCr).  The  x  and  y-axes,  respectively,  represent  the 
dilation  coefficient  m  and  the  time  delay  coefficient  n,  as  defined  in 
Eq.  (5).  The  magnitude  of  the  WT  is  color-coded  in  a  logarithmic  scale. 
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Figure  12.  Time  plot  of  for  the  awake  and  sleep  states,  respectively.  The  square  signal 
overlapped  to  the  WT,  is  high  when  the  dog  was  subjected  to  the  stimulus.  The 
WT  is  similar  in  magnitude  before  the  application  of  the  stimulus  in  both  cases.  This 
seems  to  suggest  a  background  activity  of  the  brain,  which  reacts  differenUy  to  the 
external  stimulus  depending  on  its  state  (awake  or  asleep).  In  the  case  of  awake, 
there  is  an  evident  increase  of  activity,  while  the  opposite  is  observed  in  the  sleep  state. 
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Figure  14.  A  simplified  representation  of  the  model  depicted  using  the  Mexican  hat  mother 

wavelet.  The  external  input  stimulus  generates  an  h{,t)  response  whose  amplitude  A 
is  a  function  of  the  state  of  the  subject  and  does  vary  over  time  accordingly.  The 
output  of  the  WT  signal  processor  can  be  used  to  quantify  the  time  variations  of  A(r). 
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